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LAMINAR  PIPE  FLOW  WITH  INJECTION 
AND  SUCTION  THROUGH  A POROUS  WALL* 

By 

S»  W„  Yuan*  and  A.  Finkel stein** 

Polytechnic  Institute  of  Brooklyn 

SUMMARY 

The  effect  of  injection  and  suction  at  the  wall  on  the  two-dimensional 
steady-state  laminar  flow  of  a fluid  in  a porous-wall  pipe  has  been  investigated 
in  detail  by  the  solution  of  the  Navier— Stokes  equations  in  cylindrical  coordin- 
ates. An  exact  solution  of  the  dynamic  equations,  reduced  to  a third-order  non- 
linear differential  eouation  with  appropriate  boundary  conditions,  is  obtained 0 
A perturbation  method  was  used  to  solve  the  latter  equation  for  both  small  and 
large  flows  through  the  porous  wall0 

The  velocity  components  are  expressed  as  functions  of  the  ratio  of 
velocity  through  the  porous  wall  to  the  maximum  axial  velocity  at  the  pipe 
entrance,  the  coordinates  of  the  pipe  and  the  physical  properties  of  the  fluid 0 
lhe  results  show  that  the  effect  of  injection  at  the  porous  wall  of 
the  pipe  is  to  increase  the  friction  coefficient  at  the  wall.  For  an  injection 

Q 

ratio  ^ ,01  (500  ^RQ  ^ 2500)  the  friction  coefficient  at  the  wall  is  in- 

creased by  70-05%  over  the  zero  injection  case  (Poiseuille  case). 

** 
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INTRODUCTION 


Studies  of  the  problem  of  cooling  rocket  and  jet  motors  by  the 
diffusion  of  fluids  through  porous  metal  combustion  chamber  liners  have 
been  made  by  means  of  the  investigati on  of  the  boundary  layer  behavior 
along  a porous  plate  with  fluid  injection  [1]*.  Laminar  flow  in  a two 
dimensional  channel  with  porous  walls  has  also  been  investigated  for 
extremely  small  suction  velocity  at  the  wall  [2j.  Since  the  problem 
of  the  flow  through  a porous-wall  pipe  with  injection  has  not  been 
thoroughly  investigated,  the  purpose  of  this  work  was  to  obtain  the 
basic  phenomena  of  this  type  of  flow  which  would  provide  a guidance 
for  the  investigation  of  turbulent  pipe  flow  with  injection  or  suction. 

another  important  application  of  the  results  of  this  study  is  to 
the  boundary  layer  control  for  decreasing  drag  and  increasing  lift  of 
airplane  wings.  In  this  connection  the  boundary-layer  flow  is  sucked 
through  the  surface  of  the  wing  to  a duct.  The  ensuing  flow  in  the 

duct  simulates  the  problem  of  flow  through  a pipe  with  fluid  injection 
at  the  walls . ’ 


In  the  present  investigation  an  exact  solution  of  the  Navier-Stokes 
equations  in  cylindrical  coordinates  with  injection  or  suction  as  a 
boundary  condition  at  the  wall  was  obtained.  On  the  other  hand  the 
problems  of  flow  on  a porous  flat  plate  or  curved  wall  were  previously 
made  by  the  approximate  solution  of  the  Prandtl  boundary  layer  equations. 
The  •assumptions  made  in  the  present  study  were;  (1)  the  fluid  is 
p r e s s i o 1 e , i.e.  the  mass  density  and  the  viscosity  of  the  fluid  were 
assumed  to  be  constant;  (2)  the  main  flow  was  assumed  to  be  laminar,  ancl 
the  fluid  flowing  in  the  axjal  direction  and  the  fluid  flowing  through  the 


-"'Numbers  in  brackets  refer  to  F>ibliorraohy, 
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porous  wall,  were  assumed  homogeneous;  (3)  the  maximum  axial  velocity  at 
the  entrance  of  the  porous-wall  pipe  is  equal  to  the  maximum  axial  velocity 
in  the  Poiseuille's  flow;  and  (4)  the  fluid  flowing*  through  the  porous  wall 
is  uniform  throughout. 

Fundamental  Equations 

The>  three-dimensional  steady  flow  of  a viscous  incompressible  fluid 
is  governed  by  the  following  set  of  basic,  laws  of  fluid  mechanics.  From 
Newton's  Second  law,  the  Navier-Stokes  equations  were  derived  as  follows;-. 

f = - J + U 

(1) 

From  the  principle  of  conservation  of  matter. 


v-  ? = o 


(2) 


If  a curvilinear  coordinate  system  is  introduced  with  the  origin  at 
the  center  of  the  cross-section vhere  x is  taken  in  the  direction  of  the 

flow,  r in  the  radial  direction  and  9 the  azimuthal  angle,  and  the  ordinary 

« > * 

vector  curvilinear  coordinate  transformations  are  used,  the  Navier-Stokes 
and  continuity  equations  become 


a ±2  + 

SX  ^ r 


/ -)/ 

— — /•  p , 

f d x 


. j — + y_u 

x r2  y J r x x 1 


(3) 


u -ax  dr 


I r 4 -iif  v - JC 

J + u l xx3-  Xyz  y sr  r 


uy 


3. 


£(ru)  , r)(r^)  _ 
d x + $ r ~ u 


where  u and  v represent  the  x and  r components  of  the  velocity  at  any 
point  (see  Fig.  1)  and  - 0 because  of  axially  symmetric  flow. 

The  above  equations  will  be  used  to  investigate  the  fluid  flow  in 
a circular  pipe  with  a porous  wall  through  which  uniform  fluid  injection 
or  suction  is  applied.  The  boundary  conditions  are  * 


at  r ■»  0 


' = ST  - 0 


at  r = R 


L(  — O V — - V — CcnS'f 


For  a two-dimensional  incompressible  flow  a steam  function  exists 


such  that 


Y CL  = 


- r\r  - 


and  the  continuity  equation  (5)  is  satisfied. 

c* 

*°r  a constant  fluid  injection  or  suction  at  the  porous  wall  and  the 
given  boundary  conditions,  the  following  stream  function  is  introduced; 


f = [*)■>  ov  If"/ ) 


r £ 

where  '/  « (-}  ' . The  constant  A is  determined  f r w.  the  condition  at 
7 **  ° and  x ’ 0 7,here  ihe  maximum  velocity  Ua  for  Poiseuille * s flow 
. J he  constant  B is.  determined  from  the  law  of  conservation  of 
matter.  The  stream  function  can  then  be  ux creased  as* 
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[p  _ t_-f  u0 
2 /?*) 


° + nf 


7/^ 


From  Eqs.  (5)  and  (7)  the  velocity  components  in  the  direction  of 
flow  and  the  radial  direction  are  given  by 


(9) 


t(  ~u*f  77:, + 


fit) 


V = - ^ K t 


fa) 


l/? 


(10) 


w 


U,£ 

'e  = ^ 


and  R.  = . The  function  /(/?,)  appearing  in 


here  X =»  ~^L3 
u 

tne  above  equations  is  the  only  unknown  yet  to  be  determined  in  terns  of 
the  distance  parameter  ^ Eq.  (11)  indicates  that  the  radial  velocity 
becomes  a function  of  q only.  This  is  because  of  'he  assumption  of 

constant  velocity  vQ  at  :he  wall. 

_ . «* 
ay  introducing  the  expressions  of  u and  v ‘’rom  Lqs.  (It)  and  (11)  in 

the  F/O s . v 3 ; and  (4)  there  result 


(11) 


to 


/ & = ¥ //  +4  i i]r*K(f-f(')-»(ir*n] 


(12) 


-1 

/dp 

I 


R. 


r 1 1 

I 


Since  the  right-hand  side  of  Eq.  (13)  is  a function  of 
of  both  sides  of  it  with  respect  to  x yields 


F(tj  j 

(13) 

^ only,  differention 


(14) 


0 


Hence  differentiating  Eq.  (12)  with  respect  to 


~[  (f-ff")  - “(if'+f")  ] = 0 

ciq 


which  is  to  be  satisfied  for  all  x. 


« 

Integrating  Eq*  (15),  one  obtains 


a = c 


for  \ ^ , and 


r-/r-  -ir7/^/v-A 


for  X ^ ^ ; 'where  c and  k are  the  constants  of  integration  to  be 

determined.  The  boundary  conditions  are  obtained  with  the  aid  of  Eqs.  (6), 
(10)  and  (11).  Thus: 

f(°)  = -f  ( i ) — q L<m  v'n  f (n)  — <9 

/Ui  = + 


Eq.  (16)  is  an  ordinary  non-linear  differential  equation  of  the  third 
order  which  resulted  from  the  Navier-Stokes  equations  and  the  continuity 
equation  by  the  similarity  transformation. With  the  aid  of  the  four  given 
boundary  conditions  an  exact  solution  can  be  obtained  and  the  constant  of 
integration  c determined. 


6. 


It  can  be  seen  that  the  limiting'  form  of  JSq.  (16),  by  letting  l r 
approach  to  zero,  is  the  equation  describing  a flow  through  a circular 
pipe  with  permeable  walls.  The  solution  of  this  equation  which  satisfies 
all  the  four  boundary  conditions  given  in  Eq . (18)  is  the  well-known 
Poiseuille’ s law  for  pipe  flow.  If  small  values  of  X are  treated  as 

a pertubation  parameter  a solution  of  Eq.  (16)  can  be  obtained  which  will 
be  discussed  in  the  next  section. 

On  the  other  hand,  if  iarre  values  of  A are  treated  as  a pertubation 
parameter  the  third  order  differential  equation  (17!  is  reduced  to  a second 
order  one.  The  solution  of  Eq.  (17)  can  also  be  obtained  in  the  same  manner 
since  all  four  boundary  conditions  eiven  in  Eq.  (18)  can  be  satisfied. 

SOLUTION  FTR  SKMT.T.  X 

Ihe  solution  of  Eq.  (16)  can  be  expressed  for  small  value  of  X 
by  a power  series  developed  near  A = 0 as  follows 


7. 


. "/  // 


if  + A ~ Zf  f + fo  l + fof, 


The  boundary  conditions  to  be  satisfied  by  the  f^'s  are  from  Eq.  (18) 

• •! 

f(°)  - fn(’)  = 0 ■ L,pn  Jn  f<n)  = 0 {or  an  n 

n j ^ ofn 


i<‘>  = T | 

fn  0)  =■  0 n * I ) 


The  second-order  perturbation  solution  of  Eq.  (16)  obtained  by  solving 
:-qs.  (21)  to  (23)  is  given  as  follows: 


f'V-m-i')')*  * (S-W* 


A Z 1 - 

A / 


j-4°7  ' TTz  7 - TV*?  ' 7~„7  - ^ 


- 1 - ± >■,  + — a2' 

4 .370  ' 


It  is  seen  from  the  above  equations  that  the  second-order  perturbation 
solution  is  sufficiently  accurate  even  for  \ = 1.  The  velocity  components 
ln  the^xial  and  ™4Lal  directions  are  obtained  by  substituting  Eq.  (25) 
bos.  (,10)  and  (11),  respectively  as  follows: 


8 


18  ffoo  ^ 


V .TbTc  (‘ ^ ~7Col  + 61* f-3oo n\ 


£=-■#4  0 " r?i+  +(-<■]  n-,‘- * I'n4) 


t o gee 


(•6€>J  - IS  I4*  >^rf'r)  *“ )] 


The  pressure  distribution  in  the  axial  and  radial  directions  are 
obtained  upon  the  substitution  of  Eq.  (25)  into  Eqs.  (12)  and  (13),  and 
integrating.  Then  one  obtains 


j"  (o,  o)  - -H*.  r) 


(i  + ± X - ~ \ X]  f — + 2-  A _*//_£) 

4 2?  J L //,,  p q 1 ' f?  ) 


fUj  ee  4 2;  J L y(o)  4?e  * J * 


.the  pressure  drop  in  the  flow  direction  can  bo  readily  obtained  i'roa\  £q. 


(29),  i.e. 


p(  o,  r ) - -f3  (K , ? ) 


8 r.,  3 \ //  i 1 / r < . , A x i * , 

■/Ce  4 ZJ  f (o)  *c  « ; * 


9. 


The  coeificient  of  skin  friction  at  the  wall  can  also  be  obtained  from 
Eq.  (27),  and  can  be  written 


c = __2  ~ r ‘ v 4 A.JL] f/+  ~ - — 2 1 

i ft,  Z / - -A  , JAA4  4 /?e  * '*  S4  0 J 


SOLUTION  FOR  LARGE  X 


The  solution  of  iiq.  (I?''  can  be  expressed  for  large  -valuesfof  A , 
by  a power  series  developed  near  -L  = 0 as  follows- 

A 


/ = A -1  A 4 * • ■ 


ft  — fta  4 — k i -A  kx  f • 
A 1 A 


-t  — kn 
An  n 


where  the  f^'s  and  k^'s  are  taken  to  be  independent  of  X - By 
substituting  nos.  (32)  and  (33'  into  Eq.  (17'  and  setting  all  coefficients 
of  like  powers  of  equal  to  zero,  one  obtains  the  following  set  of 

equations:- 


w - /:  - *. 


r r ,f  r * r f „ If  Jo  // 

ff  - tf of,  + f0  fx  * (7 fQ  + /0 ; = k, 


fofx  - 1 f 4 f,tL  + J,  f,  - f + (jf,  +f/  )=z  A 


I 


(35) 


10. 


The  boundary  conditions  to  be  satisfied  by  the  f *s  are  ^rom  (13) 


fn  (°)  = /nV/;  = o ^<>7  *"  a 

o ' 

/«  o)  = i | 

f,  (>)  = 0 , n z I J 

The  first-order  perturbation  of  Eq.  (17),  obtained  by  solving  the 

0 * " « 

non-linear  second-order  Eq.  (34)  and  the  linear  second  order  Eq.  (35), 
is  given  as  follows: 


J.  ^ 


where 


k - 0 (c  > Q &$  - 

A 


' ^ A J /o  -J  r 


=■  0.7  z 7 08 


JL  f S 

1 0 -/  Sin  U ^ 


/.  A / 2 


and  k - 


-1.3^53 


11. 


The  velocity  components  in  the  axial  and  radial  directions  are  obtained 
by  substituting  Eq.  (38)  into  Eqs.  (IQ)  and  (11),  respectively  as  follows: 


a 


— = F 

u,  L 


! . X x 


0-7654  4 Mi! ' 

A 


T ' 7 + tO-J 


(40) 


u ' 


2 A 

4e  JrJ 


Ci^  i'l  * -jr  fi  * ■■■  1 


(41) 


The  pressure  distribution  in  the  axial  and  radial  directions  are 
obtained  upon  the  substitution  of  Eq.  (38)  into  £qs.  (12)  and  (13),  and 
integrating 

= f C'-ceet*  +'£{]<{) 


-4  4 F -L  + % - F F / f°)J 


The  pressure 
(L2\  i.a. 


dron  in  the  flow  direction  can 


(42) 

be  readily  obtained  from  j£c. 


fuJ  ~ F~e 


l o,  6 / f.  35 


A.  /_* 

/?e  K ' K 


(43) 
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The  coefficient  of  skin  friction  at  the  wall  can  also  be  obtained  from 
Eq.  (40),  which  is 


L 0.7 8 $4 


-1 1“  <<  A 4 Iff*  -L  k,tt4 )*■■■  ] 


DISCUSSION 


The  velocity  distributions  in  the  main  flow  direction  at  an  arbitrary 
cross-section  of  the  pipe  as  calculated  from  Eq . (27)  for  A « *-1  and 
from  Eq.  (40)  for  A = 10  are  shown  in  Fig.  2.  It  was  noted  that  when 
A - 0 the  profile  becomes  Poj seuille ' s paroboloid,  and  for  A > 0 
(fluid  being  injected  through  the  wall)  the  axial  velocity  increases  and 
the  velocity  rr'nicnt  at  the  wall  increases.  For  A<  0 (fluid  being 
withdrawn  through  the  wall)  both  the  axial  velocity  and  the  velocity  gradient 
at  the  wall  decrease  as  compared  with  Foiseuille's  case.  The  above  phenomenon. 


follows  the  law  of  conservation  of  matter.  In 'the  present  case  the  radial 
velocity,  which  vanishes  in  Poiseuille's  case,  has  a finite  magnitude  except 
at  the  center  of  the  pipe  where  it  vanishes. 

In  Fig.  3 the  maximum  velocity  distributions  along  the  axis  of  the  pipe 
was  snowru  The  increase  of  maximum  velocity  with  the  increase  of  fluid 
.injection  at  the  wall  and  decrease  vith  suction  ware  illustrated  in  Fig.  4. 

It  was  interesting  to  learn  that  for  an  injection  ratio  £ -0.0/  '(Re  „ 1000) 


the  maximum  velocity  increases  about  35£  over  the  Poi seuille • 3 flov 


w case. 


13. 


The  pressure  drop  in  the  main  «*  direction  is  shown  in  Fig.  5.  It 
was  found  that  this  pressure  drop  became  appreciably  larger,  even  for  very 
email  fluid  injection  at  the  wall,  than  that  in  the  PoiseuiUcs  flow  case, 
and  it  became  appreciably  less  for  small  suction  case.  The  ratio  between 
the  pressure  drop  in  the  radial  and  the  axial  directions  is  approximately 
equal  to  the  injection  (or  suctionl  ratio  £ and  hence  the  pressure 
drop  in  the  radial  direction  can  be  neglected  in  most  practical  applications. 

One  of  the  essential  parameters  in  the  present  investigation  is  the  ' 
skin  friction  coefficient  at  the  wall.  In  Poisoulllo's  flow  the  skin 
friction  coefficient  at  the  wall,  C , has  a constant  value  of  ± 

The  wall  frictional  coefficient  as  calculated  from  Eqs.  (31)  and  (44) 
indicates  that  the  effect  of  injection  in  a pipe  flow  is  to  increase  the 
wall  frictional  coefficient  and  the  suction  to  decrease  the  wall  frictional 
coefficient.  In  a bounds ry-layer  flow  on  a porous  flat  plate  the  effect 
of  fluid  injection  at  the  wall  is  to  Increase  the  thickness  of  the  boundary 
layer  and  decrease  the  velocity  gradient  at  the  wall,  hence  the  wall  friction 
decreases  in  this  case.  On  the  other  hand  in  a pipe  flow  the  effect  of 
fluid  injection  at  the  wall  is  to  accelerate  the  main  stream  velocity  hence 

the  velocity  gradient  at  the  wall  which  determines  the  wall  friction  increases. 

For  a fluid  injection  ratio  ^ m m • _ . 

~ZT0  = •C1  the  wali  frictional  coefficient 

increases  by  a*  over  the  Poise, .ille-e  flow  case.  The  above  phenomena  were 
shown  m Figs.  6,  ? and  8.  The  comparison  of  the  variation  of  local  wall 

£ T*1  C t lOricll  COGf^T.GlPnt’  w T 4"  h f 1 i<w  j 1 • .« 

lth  fluid  ln.'Jectxon  between  the  case  of  flow  in 

3 pox' O U a 1 1,  D1  DG  and  nn  n *1  . 

'!  f-at  plat.e  was  shown  in  Fig  9, 
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FIG  2 VELOCITY  PROFILES  Vs  LENGTH  IN  RADIAL 
DIRECTION  FOR  VARIOUS  X (r«=io3,  ^-=10) 


FIG  3 MAXIMUM  AXIAL  VELOCITY  DIFFERENCE  Vs 
LENGTH  IN  THE  FLOW  DIRECTION  (Rc-IO5) 
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FIG.  5 AXIAL  PRESSURE  DROP  Vs  LENGTH 
IN  FLOW  DIRECTION  FOR  VARIOUS  X 
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FIG'?  y^lrA™K,0F  WALL  rR,CT|ONAL  COEFFICIEI 
INJECTION  AND  SUCTION  RATIO  FOR  - 

(r.=io*.  -g-=io) 
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FIG.  8 
WITH 


VARIATION  OF  WALL  FRICTIONAL  COEFFICIENT 
INJECTION  RATIO  FOR  LARGE  X (r.=ios,-^-=io) 
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VARIATION  OF  LOCAL  WALL  FRICTIONAL 
COEFFICIENT  WITH  FLUID*  INJECTION 


